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Abstract We report our recent work on mean-field potential effects on the elliptic flows of matters and antimatters 
in heavy ion collisions leading to the production of a baryon-rich matter. Within the framework of a multiphase 
transport (AMPT) model that includes both initial partonic and final hadronic interactions, we have found that 
including mean-field potentials in the hadronic phase leads to a splitting of the elliptic flows of particles and their 
antiparticles, providing thus a plausible explanation of the different elliptic flows between p and anti-p, K* and K , and 
x’ and z observed by the STAR Collaboration in the Beam Energy Scan (BES) program at the Relativistic Heavy Ion 
Collider (RHIC). Using a partonic transport model based on the Nambu-Jona-Lasinio (NJL) model, we have also 
studied the effect of scalar and vector mean fields on the elliptic flows of quarks and antiquarks in these collisions. 
Converting quarks and antiquarks at hadronization to hadrons via the quark coalescence model, we have found that the 
elliptic flow differences between particles and antiparticles also depend on the strength of the quark vector coupling in 


baryon-rich quark-gluon plasma, providing thus the possibility of extracting information on the latter's properties from 


the BES program at RHIC. 
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1 Introduction 


Heavy ion collisions at relativistic energies provide the 
possibility to study the phase structure of the strongly 
interacting matter that is described by the quantum 
chromodynamics (QCD). For the top energy available 
at RHIC and at the Large Hadron Collider (LHC), the 
produced quark-gluon plasma (QGP) is nearly baryon 
free and thus has very small baryon chemical 
potentials. According to lattice QCD calculations! ?l, 
the transition from the QGP phase to the hadronic 
matter (HM) phase in this region of the phase diagram 
is a smooth crossover without a clear phase boundary. 
The phase transition between QGP and HM is 


expected to change from the crossover to a first-order 
transition at a certain finite baryon chemical potential 
called the critical point in the QCD phase diagram”), 
To probe this region of the QCD phase diagram, the 
BES program of Au+Au collisions at much lower center 
mass energies of 7.7 GeV, 11.5 GeV and 39 GeV than 
the top energy has recently been carried out at RHIC 
by the STAR Collaboration"*!, Although no definitive 
signals for a first-order phase transition and the critical 
end point have been established, a number of 
interesting results have been observed". One of them 
is the increasing difference between the elliptic flows 
of particles and antiparticles, thus a breaking of the 
number of constituent quark number scaling of elliptic 
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flows, as the collision energy decreases. Such a 


behavior cannot be described by a simple 
hydrodynamic or hadronic cascade model even if the 
quark coalescence is considered during the 
hadronization of produced QGP"! Several theoretical 
attempts have been made to explain this surprising and 
interesting experimental result. These include the 
different elliptic flows of transported and produced 
partons during the initial stage of heavy ion 
collisions''"!; the chiral magnetic effect induced by the 
strong magnetic field in non-central heavy ion 


[12] the hadronic"! and the partonic!“ 


collisions 
mean-field effects, and the local thermal and baryon 
chemical equilibrium effect!?!. In the present talk, we 
review the results on the mean-field effects on the 
elliptic flows of particles and antiparticles in heavy ion 


collisions at BES energies. 


2  Hadronic 
elliptic flows 


mean-field potentials and 


In Ref.[13], we have studied quantitatively the elliptic 
flows of particles and their antiparticles at BES 
energies by extending the AMPT modell?! to include 
the potentials of baryons, kaons, and pions as well as 
their antiparticles. The AMPT model is a hybrid model 
with the initial particle distributions generated by the 
(HIJING) 
via the Lund string fragmentation model. In 


Heavy-Ion Jet Interaction Generator 


model"! 
the string melting version of the AMPT model, which 
is used in Ref.[13], hadrons produced from excited 
strings in the HIJING model are converted to their 
valence quarks and antiquarks, and their evolution in 
time and space is then modeled by Zhang's parton 
cascade (ZPC) model'?, Different from previous 
applications of the AMPT model for heavy ion 
collisions at higher energies, the parton scattering 
cross section and the ending time of the partonic stage 
are adjusted in Ref.[13] to approximately reproduce 
measured elliptic flows and the hadronic energy 
density (—0.30—0.35 GeV/fm?) at the extracted baryon 
chemical potential and temperature at chemical freeze 


out using the statistical model!'?l, 


Specifically, the 
parton scattering cross section is taken to be isotropic 
with the values 3 mb, 6 mb and 10 mb and the ending 


times of the partonic stage to be 3.5 fm/c, 2.6 fm/c and 


2.9 fm/c for collisions at the center of mass energies of 
7.7 GeV, 11.5 GeV and 39 GeV, respectively. At 
hadronization, quarks and antiquarks in the AMPT 
model are converted to hadrons via a spatial 
coalescence model, and the scatterings between 
hadrons in the hadronic stage are described by a 
relativistic transport (ART) model”! that has been 
extended to also 


include — particle-antiparticle 


annihilations and their inverse reactions. 


2.  Hadronic potentials 


It is known from heavy ion collisions at lower 
collision energies at SIS/GSI and AGS/BNL that the 
elliptic flow of nucleons is affected not only by their 
scattering but also by their mean-field potentials in the 


hadronic matter?! 


. This is because particles with 
attractive potentials are more likely to be trapped in 
the system and move in the direction perpendicular to 
the participant plane while those with repulsive 
potentials are more likely to leave the system and 
move along the participant plane, thus reducing and 
enhancing their respective elliptic flows. Also, the 
potentials of a particle and its antiparticle are different, 
and they generally have opposite signs at high 
densities"?! As a result, particles and antiparticles 
are expected to have different elliptic flows in heavy 
ion collisions when the produced matter has a nonzero 
baryon chemical potential. Furthermore, the difference 
between the potentials of a particle and its antiparticle 
diminishes with decreasing baryon chemical potential, 
so their elliptic flows are expected to become similar 
in higher energy collisions when more antiparticles are 
produced. These effects are all consistent with what 
have been seen in the experimental data from the BES 
program. 

For the nucleon and antinucleon potentials, 
we take them from the relativistic mean-field model 
used in the Relativistic Vlasov-Uehling-Uhlenbeck 
transport model*! in terms of the nucleon scalar and 
vector self-energies in a hadronic matter. Because of 
the G-parity, nucleons and antinucleons contribute 
both positively to the scalar self-energy but positively 
and negatively to the vector self-energy, respectively. 
Since only the light quarks in baryons and antibaryons 
contribute to the scalar and vector self-energies in the 
mean-field approach, the potentials of strange baryons 
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and antibaryons are reduced relative to those of 
nucleons and antinucleons according to the ratios of 
their light quark numbers. For the kaon and antikaon 
potentials in nuclear medium, they are also taken from 
Refs.[24,25] based on the chiral effective Lagrangian 
that fits the 
antikaon-nucleus scattering. The resulting potential 


empirical data on kaon- and 
is then repulsive for a kaon and attractive for an 
antikaon. The pion potentials are related to their 
self-energies and have been calculated in Ref.[26] 
from the pion-nucleon s-wave interaction up to the 
two-loop order in chiral perturbation theory. In 
asymmetric nuclear matter, this leads to a splitting of 
the mean-field potentials for positively and negatively 
charged pions. 

In the absence of antibaryons, the nucleon 
potential is slightly attractive while that of antinucleon 
is strongly attractive, with values of about —60 MeV 
and —260 MeV, respectively, at normal nuclear matter 
density po=0.16 fm^. The latter is similar to that 
determined from the non-linear derivative model for 


] 


. S ; 27 ; 
small antinucleon kinetic energies! and is also 


consistent with those extrapolated from experimental 


data on antiproton atoms 931 


In antibaryon-free 
matter, the K^ potential is slightly repulsive while the 
K potential is deeply attractive, and their values at po 
are about 20 MeV and —120 MeV, respectively, similar 


34-37] 


to those extracted from experimental data! and 


used in the previous study”®. 


For pions in 
neutron-rich nuclear matter, the potential is weakly 
repulsive and attractive for z and T, respectively, and 
the strength at po and isospin asymmetry 
Ô=(Pn-Pp)/(PrtPp)=0.2 is about 14 MeV for z and 


—] MeV for z 9l. 


22 Particle and antiparticle elliptic flows 


The differential elliptic flows of p, K', and z as well 
as their antiparticles with respect to the participant 
plane with and without hadronic potentials at three 
different BES energies from the string melting AMPT 


model are shown in Fig.1. Without hadronic potentials 
elliptic flows from the AMPT model are similar for 
particles and their antiparticles, including hadronic 
potentials which slightly increases the p and anti-p 
elliptic flows at p; «0.5 GeV/c, while reduces slightly 
(strongly) the p (anti-p) elliptic flow at higher pr. 
Hadronic potentials also increase slightly the elliptic 
flows of K` while mostly reduce that of K`. In addition, 
the effect from potentials on the elliptic flow decreases 
with increasing collision energy, which is consistent 
with the decreasing net baryon density of produced 
hadronic matter with increasing collision energy. The 
difference between the differential elliptic flows of p 
and anti-p, and between those of K^ and K below 
collision energy of 11.5 GeV are qualitatively 
consistent with experimental data"! while that of x 
and z is small in all three energies due to the small 
isospin asymmetries of produced matter. 

Our results for the relative printegrated v2 
difference between particles and their antiparticles, 
defined by [v2(P)-v2(anti-P)]/v2(P), with and without 
hadronic potentials are shown in Fig.2. These 
differences are very small in the absence of hadronic 
potentials, including hadronic potentials increases the 
relative v; difference between p and anti-p and 
between K' and K` up to about 30% at 7.7 GeV and 
20% at 11.5 GeV but negligibly at 39 GeV. These 
results are qualitatively consistent with the measured 
values of about 63% and 13% at 7.7 GeV, 44% and 3% 
at 11.5 GeV, and 12% and 1% for the relative v2 
difference between p and anti-p and between K^ and 
K, respectively. Similar to experimental data, the 
relative elliptic flow difference between z' and z is 
negative at all energies after including their potentials, 
although ours have smaller magnitudes. We have also 
found that, as seen in experiments", the relative v? 
difference between A hyperon and anti-A is smaller 
than that between p and anti-p, because the A (anti-A) 
potential is only 2/3 of the p (anti-p) potential. 
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Fig.1 (Color online) Differential elliptic flows of mid-rapidity (|y|<1) p and anti-p [(a), (b), (c)], K" and K` [(d)(eYf)], and z^ and x" 
[(g)(h))] with and without hadronic potentials in Au+Au collisions at b = 8 fm and center of mass energy of 7.7 GeV(a)(d)(g), 11.5 
GeV(b)(e)(h), and 39 GeV (c)(f)(i) from the string melting AMPT model. The solid (dashed) curves represent the results of particles 
(anti-particles) while the thick (thin) curves represent the results with (without) potentials. 
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Fig.2 (Color online) Relative elliptic flow difference between 
mid-rapidity (|y|«1) p and anti-p, K' and K , and z' and z with 
(open symbols) and without (solid symbols) hadronic potentials 
in Aut Au collisions at 5-8 fm and the center of mass energy of 
7.7 GeV, 11.5 GeV and 39 GeV from the string melting AMPT 
model. 


3  Partonic potentials and elliptic flows 


3.1 Nambu-Jona-Lasinio model 


To study the effect of partonic mean fields on parton 


elliptic flows, we have used in Ref[14] the NJL 


39,40] 


modelP?^" particularly the one for three quark 


flavors". Like particles and antipartcles in baryon-rich 
hadronic matter, quarks are affected by attractive 
scalar and repulsive vector fields, while antiquarks are 
affected by attractive scalar but attractive vector fields 
in baryon-rich quark matter. The value of the attractive 
scalar mean field is related to the difference between 
the current and the constituent quark masses, which 
has a maximum value of about —300 MeV. The 
magnitude of vector mean field depends on the 
product of the quark vector coupling and the 
net-baryon density. For the vector coupling gj&17 
MeV fim? obtained from the Fierz transformation of the 
quark scalar interaction in the NJL Lagrangian and 
used in calculating the results shown in Figs.3 and 4 of 
Section 3, the vector mean field has a magnitude of 
=17 MeV if the net-baryon density is 1 fm^, leading 
to a difference of «34 MeV in the quark and antiquark 
vector mean fields. 

Similar to that for the hadronic matter ^ ^? 
the time evolution of the partonic matter produced in 
relativistic heavy ion collisions can be described by 
the transport equation for the parton phase-space 
distribution function. 
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3.2 Quark and antiquark elliptic flows 


For the initial quark and antiquark rapidity and 
in AutAu 
collisions at center of mass energy of 7.7 GeV and at 


transverse momentum distributions 
impact parameter b=8 fm, we use the valence quarks 
and antiquarks converted from hadrons that are 
obtained from the HIJNIJ modell?! through the Lund 
string fragmentation as implemented in the AMPT 


model with string melting! 9. 


Including partonic 
mean fields and using a constant isotropic parton 
Scattering cross section of 5 mb, we have solved the 
partonic transport equation with the test particle 
method’. Figure 3 shows the transverse momentum 
dependence of quark and antiquark elliptic flows at the 
end of the partonic phase, which is about 1.9 fm/c after 
the start of the partonic evolution when the energy 
density in the center of produced quark matter 
decreases to 0.8 GeV/fm?. Shown are the four cases of 
mean field, scalar and 


including only scalar 


time-component vector mean fields, scalar and 
space-component vector mean fields, and scalar and 
full vector mean fields. For the scalar mean field, 
which is attractive for both quarks and antiquarks, it 
leads to a similar reduction of the quark and antiquark 
elliptic flow as first found in Ref.[45]. The vector 
mean field, on the other hand, has very different 
effects on quarks and antiquarks in the baryon-rich 
QGP as it is repulsive for quarks and attractive for 
antiquarks. The time component of vector mean field 
turns out to have the strongest effect, resulting in a 
significant splitting of the quark and antiquark elliptic 
flows as a result of enhanced quark elliptic flow and 
flow. The 
component of vector mean field has, however, an 


suppressed antiquark elliptic space 
opposite effect; it suppresses the elliptic flow of 
quarks and enhances that of antiquarks, although 
relatively small and appearing later in the partonic 
stage compared to that of the time component of 
vector mean field. 

For the case including the scalar and both 
components of vector mean field, the difference 
between the integrated elliptic flow of up and down 
quarks and their antiquarks is about 60% of the 
integrated elliptic flow of up and down quarks, and 
that between strange quarks and anti-strange quarks is 


about 2996 of the integrated elliptic flow of strange 
quarks. 
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Fig.3 (Color online)Elliptic flows of light and strange quarks 
and antiquarks at midrapidity (|v|«1) as functions of transverse 
momentum at hadronization for the cases of including only 
scalar mean field (S) (a), scalar and time-component vector 
mean fields (S+Vo) (b), scalar and space-component vector 
mean fields (S+V;) (c), and scalar and full vector mean fields 
(S*Vg*V;) (d). 


3.3 . Hadronic elliptic flows 


To study how different quark and antiquark elliptic 
flows are reflected in those of produced hadrons, we 
have used the coalescence model to convert them to 
46] In this model, the 
probability for a quark and an antiquark to form a 


hadrons at hadronization 


meson is proportional to the quark Wigner function of 
the meson with the proportional constant given by the 
statistical factor for colored spin-1/2 quark and 
[47,48] 
meson .The 


probability for three quarks or antiquarks to coalesce 


antiquark to form a colorless 
to a baryon or an anti-baryon is similarly proportional 
to the quark Wigner function of the baryon. Since x 
and z elliptic flows are affected similarly by the 
partonic mean fields as a result of similar elliptic flows 
for u and d quarks and their antiquarks, we have 
considered only p, K, and A and their antiparticles. 

In Fig.4, we show by solid and dashed lines, 
respectively, the elliptic flows of p and anti-p (left 
panel), A and anti-A (middle panel), and K^ and K` 
(right panel), at hadronization as functions of 
transverse momentum. It is seen that the quark 
coalescence leads to a larger hadron elliptic flow than 
the quark elliptic flow at same transverse momentum. 


Furthermore, the elliptic flows of p, A, and K are 
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respectively larger than those of anti-p, anti-A, and KË, 
leading to the relative differences between their 
integrated elliptic flows of about 45, 40, and —6.0%, 
respectively, as shown by solid symbols in Fig.5 for 
the quark vector coupling gy = 17 MeV fm’, compared 
with 63414, 54427, and 1342% 
experiments shown by open symbols in the left side of 
Fig.5. 


measured in 


0.15 


n 


Pr | GeV'c? 


Fig.4 (Color online) Elliptic flows of mid-rapidity (|v|«1) p 
and anti-p (left panel), A and anti-A (middle panel), and K^ and 
K (right panel), at hadronization as functions of transverse 
momentum. 


100 


80 A-A (Exp.) 


60 | P-P (Exp.) 


“| 


(vn) v. | % 


0.0 0.2 0.4 0.6 0.8 
gy G 


Fig.5 (Color online) Relative difference between integrated 
elliptic flow of mid-rapidity (|v|«1) p and anti-p (solid squares), 
A and anti-A (solid circles), and K* and K (solid triangles) at 
hadronization for several values of quark isoscalar vector 
coupling. Experimental data from Refs.[8,9] are shown by open 
symbols in the left side. 

The dependence of the relative difference 
between integrated particle and antiparticle elliptic 
flows on the vector coupling gy is shown in Fig.5. 
Besides the value g;)/G=0.33, corresponding to gy-17 
MeV fm’, two other values of 0.165 and 0.73?! are 


also used. The relative integrated elliptic flow 


differences are seen to increase almost linearly with 


the strength of the vector coupling. 


4 Conclusion 


We have studied the elliptic flows of p, K’, x and 
their antiparticles in heavy ion collisions at BES 
energies by extending the string melting AMPT model 
to include their mean-field potentials in the hadronic 
stage!!! Because of the more attractive anti-p than p 
potentials, the attractive K and repulsive K^ potentials, 
and the slightly attractive z' and repulsive az 
potentials in the baryon- and neutron-rich matter 
formed in these collisions, smaller elliptic flows are 
obtained for anti-p, K , and a than for P, K', and 7. 
Also, the differences between the elliptic flows of 
particles and their antiparticles are found to decrease 
with increasing collision energy as a result of 
decreasing baryon chemical potential of produced 
hadronic matter. Although our results are qualitatively 
consistent with the experimental observations, they 
somewhat underestimate the relative elliptic flow 
difference between p and anti-p as well as that 
between z and z and overestimate that between K^ 
and K. 

We have also studied the effect of partonic 
mean fields on the elliptic flows of quarks and 
antiquarks in a baryon-rich quark matter by using a 
transport model based on the NJL modell'“, Although 
the scalar mean field, which is attractive for both 
quarks and antiquarks, leads to a similar reduction of 
the quark and antiquark elliptic flow, the vector mean 
field, on the other hand, results in a significant 
splitting of the quark and antiquark elliptic flows as a 
result of enhanced quark elliptic flow and suppressed 
antiquark elliptic flow. Using the quark coalescence 
model, we have further studied the elliptic flows of p, 
A, and K^ and their antiparticles produced from the 
and found that the 
differences between particle and antiparticle elliptic 


baryon-rich quark matter 


flows are appreciable as a result of the different quark 
and antiquark elliptic flows. The magnitude of the 
relative integrated elliptic flow difference between 
particles and their antiparticles depends on the strength 
of the vector coupling. For the small vector coupling 
obtained from the Fierz transformation of the scalar 
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interaction in the NJL model (g,/G=0.33), the relative 
elliptic flow difference due to partonic mean fields is 
about one third larger for p and anti-p, close to a factor 
of two larger for A and anti-A, and about a factor of 
five smaller but with an opposite sign for K” and K, 
compared to corresponding values due to the effect of 
hadronic potentials given in Ref.[13]. 

As to the elliptic flows of z' and 7, they are 
the same in our study reported in Ref.[14] as we have 
not included the isovector part of partonic mean fields. 
As shown in our study reported in Ref.[13], including 
the isovector hadronic mean fields indeed leads to a 
splitting of the elliptic flow of z and z, although a 
factor of five smaller than the measured value. It is 
thus of great interest to study the effect of quark 
isovector mean field on the z^ and z difference and 
compare it with that due to the chiral magnetic effect 
suggested in Ref.[12]. 

If the effects of partonic and hadronic 
mean-field potentials were additive, which certainly 
needs to be checked in future studies, then the final 
relative elliptic flow differences between p and anti-p, 
A and anti-A, and K' and K would be close to those 
measured in the BES experiments at RHIC if a small 
quark vector coupling is used. The latter is consistent 
with that extracted in Refs.[50-52] from the baryon 
number susceptibilities measured in lattice QCD. A 
more quantitative study that includes both partonic and 
hadronic mean-field effects as well as the chemical 
reactions in the partonic phase thus provides the 
possibility of extracting information on the strength of 
quark vector interaction in baryon-rich QGP and thus 
its equation of state from the BES program at RHIC. 
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